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We investigate a model of modified gravity recovering the modified 
Newtonian dynamics (MOND) in the non-relativistic limit, based on the 
introduction of a preferred time foliation violating Lorentz invariance in 
the weak-field regime. Lorentz-invariance violation has been studied in the 
framework of Einstein-asther theory, the generalization of which, known 
as non-canonical Einstein-aether theory, having been proposed as a rel- 
ativistic formulation of MOND. Our model can be seen as a minimal 
specialization to the hypersurface orthogonal case, which allows a differ- 
ent interpretation in terms of the preferred time : it can be either treated 
as a dynamical scalar field in a 4D formulation, or chosen as the time 
coordinate in a 3+1 formulation. We discuss the equivalence of the two 
points of view and the non-relativistic limit of the model. 

1 Motivation 

As an alternative to the standard A-CDM model, performing well in cosmology 
and on large scales but leaving open questions at galactic scales, the Modified 
Newtonian Dynamics (MOND) [12] intends to answer the missing mass problem 
by a phenomenological modification of gravity in the weak-field, large distances 
regime. It has proven successfull concerning galactic observations, while needing 
some extension at the scale of clusters of galaxies and for cosmology (for reviews, 
see P31[S]). We will adopt for MOND its modified Poisson-equation version [2J: 



where p is the ordinary matter density, (f> is the Newtonian potential, g = V</> is 
the gravitational field, ao = 1.210 -10 ms~ 2 is a universal constant acceleration 
scale, and /i is an interpolation fonction of the ratio x = |flf|/ffl0) with asymptotic 
behaviour : n(x) — > 1 [x — > +oo] and fJ,(x) ~ x [x — > 0]. 
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This modified Poisson equation is a purely non-relativistic formula, while a 
relativistic theory is needed to address issues concerning cosmology and gravi- 
tational lensing. Several relativistic extensions for MOND have been proposed 
in the past, including a tensor- vector-scalar theory (TeVeS) [HUE], a bimetric 
theory [T3] , non-canonical Einstein- iEther theories [TH] [5] , and a modified dark 
matter theory [3J. 

With a completely different motivation, the Hofava-Lifshitz approach [9] in- 
troduced a Lorentz invariance breaking preferred foliation of spacetime to build 
a power-counting renormalizable theory of gravity in the strong-field regime. 
This approach has been improved and extended in [5]. 

In this contribution, which is a summary of our recent paper [3] (see also 
[15] ) to which we refer for more details, we investigate a specific relativistic 
extension of MOND, that can be either understood as a Lorentz- invariance 
breaking theory using a preferred time foliation, or as a hypersurface-orthogonal 
restriction of non-canonical Einstein- iEther theories. 



2 Non-canonical Einstein-iEther theories 

Einstein- iEther theories have been introduced as a phenomenological approach 
to Lorentz-invariance violation |11] . A preferred direction of time is described 
by a dynamical, timelike unit vector n^. The corresponding Lagrangian density 
is usually written as0 



C 



M 
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R + K + X^n^ + l) 



(2) 



where R is the Ricci scalar, A is a Lagrange multiplier enforcing the normaliza- 



T, and where K. represents the most general Lagrangian 



tion condition n^n^ 

density that is quadratic in the derivatives of n 

K = Kf a,p "V ll n l ,V v n a , (3a) 
K^"° = c l9 ^g»° + c 2 g^g ua + c 3 g^g^ + c A n»n v g"° , (3b) 

where c\, C2, C3 and C4 are dimensionless constants, left unspecified at this stage. 
Constraints in the Solar System computed from Parametrized Post-Newtonian 
parameters have been obtained in [7]. 

The non-canonical generalization of the previous action has been proposed 
in [18] |8], as a relativistic extension of MOND. It consists of replacing the aether 
action JC by a free function F(K,), i.e. of considering the Lagrangian 



-'non-canonical JE 
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R + F(JC) + A(n"n^ + 1) 



(4) 



More generally, one could imagine introducing several arbitrary functions F\, 
F2, F3 and F4 corresponding to the various terms in ((3J. In the non-relativistic 



1 We use the ( h H — h) convention for the metric, and often set G •■ 
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limit, specifying the behaviour of the fonction F(JC) around JC = allows to 
recover the modified Poisson equation JT]). The class of cosmologies that are 
possibly produced by this type of models has been extensively studied in [15] . 



3 From the iEther to the Khronon 

Instead of directly using a timelike dynamical vector field in the action, 
Lorentz invariance violation can be realized by introducing a preferred foliation 
of spacetime. The link between these two formulations has been studied in 
[TU] . The foliation consists of hypersurfaces of constant value of a scalar field r, 
which we will call the "Khronon" field following [5]. It defines a hypersurface- 
orthogonal unit vector field rip, whose relation to the Khronon field r is: 

Up = -NdpT , with N = - 1 . (5) 

^/-gP <T d p Td (T T 

We see that the condition for the Khronon field to define a well-behaved space- 
time foliation is that its gradient must remain timelike everywhere. Recall that, 
in general, a unit timelike aether vector does not define a foliation of spacetime; 
the condition for it to be hypersurface-orthogonal is n[p\7 u n p ] = (the so-called 
Frobenius theorem). Notice also that the fundamental ingredient of the theory 
is now the scalar field r, while in Einstein- iEther theory it is the aether vector 
field, with three degrees of freedom (one being suppressed by the unit-norm 
constraint) instead of one. An advantage of the formulation ([5]) is that the nor- 
malization condition is automatically satisfied and we do not need a Lagrange 
multiplier in the action. 

Given this spacetime foliation, further geometrical definitions are: 

7// = + n ^ n " > K iw = p n v > a n = n v V u np , (6) 

where 7^" is the projector orthogonal to n M , Kp V is the hypersurface's extrinsic 
curvature tensor and is the spacelike four-acceleration of the congruence 
orthogonal to the hypersurface with velocity n M . We also define the projected 
covariant derivative operator as, for instance: DpV u = 7 /J p 7 (T 1 ' V pV 1 . Those 
definitions could be done formally in the non-hypersurface-orthogonal case, but 
the following important relation is specific to this case: 

a fl = D fl lnN (7) 

and, in addition, the extrinsic curvature is symmetric: — K u ^. 

We are now able to rewrite the four terms written in Eq. <j3j) as the basic 
ingredients of the Lagrangian density of the Khronon field, expressed in terms 
of extrinsic curvature and acceleration: 

Term c x : V ^V^n" = K^K^ - a 2 , (8a) 

Term c 2 : (V^) 2 = K 2 , (8b) 

Term c 3 : V.^V/ = K^K^ , (8c) 

Term c 4 : n^rfV ^n p V v n p = a 2 , (8d) 
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where a 2 = a^a^ and K = KJ* is the trace of the extrinsic curvature. We see, 
as was pointed out in [10] , that only three out of these four terms are actually 
independent in the hypersurface-orthogonal case. We can go further by inves- 
tigating the fates of the extrinsic curvature and acceleration in two regimes of 
interest, the non-relativistic or post-Newtonian limit and the homogeneous and 
isotropic cosmology. In the post-Newtonian limit and in adapted coordinates, 
i.e. when the time coordinate t is identified with t, one obtains (see below): 

a, = ^d t cj) + 0(4) , and Kij, K = 0(3) , (9) 

with (f> the Newtonian potential and where 0(n) means terms of order at least 
l/c™. On the other hand, when considering a perfectly homogeneous and 
isotropic (FLRW) Universe, where the foliation is identical to the cosmic time 
foliation, we have: 

a ll = , and KijK ij = 3H 2 , K 2 = 9H 2 , (10) 

where H is the standard Hubble parameter. 

This means that the c\ and C4 terms are the only ones to contribute in the 
non-relativistic regime, and that the C4 term vanishes in the usual cosmological 
background. However, notice that the c\ term switches sign between the two 
regimes, whereas the C4 term is positive in any case. In the model discussed 
below, we use only this a 2 term as the basic ingredient to recover MOND in 
the non-relativistic limit as an alternative to dark matter at the galactic scale, 
leaving the cosmology aside. 



4 Specific example of MOND theory 
4.1 Covariant formulation 

In the 4D covariant point of view, the model constitutes a modification of 
General Relativity (GR) by the introduction of an additional scalar field, the 
Khronon r. We introduce in the action a free function of the norm of the 
acceleration a: 

C = ^-[R-2f(a)]+C m [g^^}. (11) 

This choice corresponds to the C4 term in the Einstein- iEther action. We also 
assume the standard coupling of matter fields to the metric. The function 
f(a) is to be specified later, in order to recover MOND in the weak-field limit 
and GR in the strong-field regime. Variation of the action with respect to the 
metric and to the r field yields the field equations: 

G^ v + f(a)g^ + 2nW p [x(a)a p ] - 2 X {a)a»a u = 8tt , (12a) 

V M KV„( X (aK)] - \n u V v f + X (o) aWK^ , (12b) 
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with T^" the matter stress-energy tensor, and where we used the short- hand 
notation x( a ) = f'( a )/^ a - The second equation, which we call the T-equation, 
is of fourth order in the derivatives. We will see however that, in adapted 
coordinates, it becomes of first order only in time derivatives of geometrical 
quantities. If, from the first equation, which we call the modified Einstein 
equation, we define an equivalent stress-energy tensor T£ v for the Khronon 
field, one can see that the T-equation is in fact equivalent to V^T^ = 0. Hence, 
because of the Bianchi identity, the matter conservation equation V„T Miy = 
and the modified Einstein equation together contain the r-equation. 



4.2 3+1 formulation 

We now write the 3+1 formulation of the theory in adapted coordinates, choos- 
ing t = r. With standard definitions for the lapse N, the shift Ni and the spatial 
metric 7-y, the 3+1 parametrization of the metric reads 



ds 2 = -{N 2 - NiN^dt 2 + 2N i dtdx i + 'y ij dx i dx j 



(13) 



In these adapted coordinates, we have n M = {—N, 0) and a% = di In TV, with ./V 
being now a geometrical quantity. The 3+1 form of the action is then: 

The r field is not a dynamical field anymore — it has been absorbed into the time 
coordinate. Its contribution to the action is purely geometric, f(a) depending 
now only on the gradient of the lapse N. This term explicitly breaks Lorentz 
invariance. Varying the action with respect to N, Ni and 7^ , we obtain: 



TZ + K' 



16ne , (15a) 
D s (K ij - Y j K) = -8nJ l , (15b) 



+2K ik K\ - KK ij - -Y J (K kl K, 



^ k K) 



K 



-— (D'D^N - f j D k D k N) - 2xa i a j + f'f- 
where S£ m /SN = -y/js, SC m /5Ni = and 5C m /8-f rj 



-- SttV 3 , (15c) 
N^T ij /2. 



4.3 Equivalence between the formulations 

The two formulations of our model seem quite different, since the r field is a 
dynamical field in the 4D point of view while it becomes a mere time coordinate 
in the 3+1 point of view. It has been pointed out in |10] that the fact that the r- 
equation is contained in the modified Einstein equation and the conservation of 
the matter stress-energy tensor allows one to consistently treat it as a coordinate. 
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We explicitly checked this equivalence between the two formulations at the 
level of the field equations, by verifying that the 3+1 projection (after variation) 
of the 4D equations was in agreement with the 3+1 ones, obtained by the vari- 
ation with respect to geometrical quantities after projection. We also checked 
that the 3+1 projection of the r-equation, namely 

+ NK (Afro*) + Xo 2 ) - A \ -/'«' A',., = , (16) 

where D t = d t - N k D k , was indeed contained in the field equations (1151) and 
the 3+1 writing of the conservation of matter stress-energy tensor. Notice that 
this equation (|16|) has now become of first order only in time derivatives of 
geometrical quantities. 



A 



Xa - 77 



4.4 Recovering MOND in the non-relativistic regime 

For a system at rest with respect to the preferred frame, we may write, in the 
non-relativistic (NR) limit, restoring the c factors: 

7V = l + ^+0(4), JVi = 0(3), 7tf=*tf (l-|r)+0(4). 

Here <j> and tp are the usual Newtonian potentials, and we recall that <Zj = 
di4>/c 2 + 0(4). Combining the 3+1 field equations (|15[) . we first obtain the 
equality of the Newtonian potentials, <j> = ip + 0(2), crucial for the dark matter 
seen by gravitational lensing |14j . and we get a modified Poisson-like equation, 

A [(1 + %K] + / + a 2 - ±D t K - K^Kij =Ait(e+ ^N t J* + T) , (18) 

which reduces in the NR limit to V • [(1 + \) V0] = AnGp + (2), with p 
the ordinary rest-mass density of matter. We see that there is a one-to-one 
correspondence between the MOND p function and the function / in the action, 
namely p = 1 + x where x( a ) = /'(°)/2a. Requirements on / to recover MOND 
are as follows (Aoo and Ao being two constants): 

f(a) ~ Aqo in the strong-field regime a — > oo, (19a) 
2a 3 /a 4 \ 

/(a) = A - a 2 H h I in the weak-field regime a -> 0. (19b) 



3ao V a 



Note that the usual strong- field condition would be f'(a)/a — > (so that we re- 
cover the standard Poisson equation), but here we imposed a stronger condition 
when a — > oo in order to recover exactly GR with a cosmological constant Aqo 
in the strong-field regime. 



4.5 Solar System effects 

We can get a rough estimate of the smallness of MOND effects in the Solar 
System, assuming that (|19aj) holds, which as we said is in fact more restrictive 



G 



than the usual MOND requirement x( a ) ~ * 0. To this end we write (with c = 1) 



/(a) ~ Aoo + ka\ 




(when a — » oo) , 



(20) 



with fc a dimensionless number of order one. This constant is fixed in order of 
magnitude by the zero-point of the function / and the numerical coincidence 
A ~ ajj. Then, translating the MOND equation in spherical symmetry as (1 + 
x)g = ,9n, with gN the Newtonian acceleration, and inverting, one obtains 



Defining the MOND transition radius Rq by ao = GM/Rq, we have that 
(a /.g N ) 2+Q = {r/R a ) i+2a . Since numerically R ~ 7100 AU, we see that within 
a Neptune orbit (~ 30 AU), the relative effect is at most of the order of 10~ 12 
for a close to 0, and is at most 10~ 15 for a typical value a ~ 1. Relativistic 
corrections in the Solar System are of the order of (v/c) 2 ~ 10~ 8 , and typical 
constraints on the values of PPN parameters |17) correspond to deviations from 
GR of the order of 10 -8 x 10~ 4 = 10~ 12 . Our requirement l|20p for recovering 
GR implies that the theory passes the Solar-System tests. 

5 Conclusion 

We investigated a simple model for a relativistic extension of MOND based 
on a preferred time foliation. We discussed how it relates to the more general 
framework of non-canonical Einstein- /Ether theories and why it selects out of it 
the minimal ingredients required for recovering MOND in the non-relativistic 
limit. We analyzed the two points of view for the model, 4D or 3D, and discussed 
their equivalence. We investigated the non-relativistic limit and estimated the 
size of Solar-System effects. Further work is needed to study the cosmology of 
the model, and possibly to extend it in order to account for the missing mass 
problem at large cosmological scales. 
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